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1 Introduction and Historical Perspective 



The theory of nuclear forces has a long history (cf. Table [TJ. Based upon 
the seminal idea by Yukawa [I], first field-theoretic attempts to derive the 
nucleon-nucleon (NN) interaction focused on pion-exchange. While the one- 
pion exchange turned out to be very useful in explaining NN scattering data 
and the properties of the deuteron [2j , multi-pion exchange was beset with 
serious ambiguities Oil]. Thus, the "pion theories" of the 1950s are gen- 
erally judged as failures — for reasons we understand today: pion dynamics 
is constrained by chiral symmetry, a crucial point that was unknown in the 
1950s. 

Historically, the experimental discovery of heavy mesons [S] in the early 
1960s saved the situation. The one-boson-exchange (OBE) model [6] [7] 
emerged which is still the most economical and quantitative phenomenol- 
ogy for describing the nuclear force j8j [9]. The weak point of this model, 
however, is the scalar-isoscalar "sigma" or "epsilon" boson, for which the 
empirical evidence remains controversial. Since this boson is associated with 
the correlated (or resonant) exchange of two pions, a vast theoretical effort 
that occupied more than a decade was launched to derive the 2-7r-exchange 
contribution to the nuclear force, which creates the intermediate range at- 
traction. For this, dispersion theory as well as field theory were invoked 
producing the Stony Brook [10], Paris [H][l2], and Bonn [3 [13] potentials. 

The nuclear force problem appeared to be solved; however, with the 
discovery of quantum chromo dynamics (QCD), all "meson theories" were 
relegated to models and the attempts to derive the nuclear force started all 
over again. 

The problem with a derivation from QCD is that this theory is non- 
perturbative in the low-energy regime characteristic of nuclear physics, which 
makes direct solutions impossible. Therefore, during the first round of new 
attempts, QCD-inspired quark models [13] became popular. These models 
are able to reproduce qualitatively and, in some cases, semi-quantitatively 
the gross features of the nuclear force |15| 116) . However, on a critical note, 
it has been pointed out that these quark-based approaches are nothing but 
another set of models and, thus, do not represent any fundamental progress. 
Equally well, one may then stay with the simpler and much more quantita- 
tive meson models. 

A major breakthrough occurred when the concept of an effective field 
theory (EFT) was introduced and applied to low-energy QCD. As outlined 
by Weinberg in a seminal paper [T7] , one has to write down the most general 
Lagrangian consistent with the assumed symmetry principles, particularly 



3 



Table 1: Seven Decades of Struggle: The Theory of Nuclear Forces 



1 Q'i'^ 

±yoo 


x iiKciwM. ivieson ± neory 




T/ie Pion Theories 


iyou s 


One-Pion Exchange: o.k. 




Multi-Pion Exchange: disaster 




Many pions = multi-pion resonances: 


±yt>u s 


er, p, a;, ... 




The One-Boson-Exchange Model: success 




Refined meson models, including 


1970's 


sophisticated 2tt exchange contributions 




(Stony Brook, Paris, Bonn) 




Nuclear physicists discover 


1980's 


QCD 




Quark Cluster Models 




Nuclear physicists discover EFT 


1990's 


Weinberg, van Kolck 


and beyond 


Back to Pion Theory! 




But, constrained by Chiral Symmetry: success 



the (broken) chiral symmetry of QCD. At low energy, the effective degrees 
of freedom are pions and nucleons rather than quarks and gluons; heavy 
mesons and nucleon resonances are "integrated out" . So, the circle of his- 
tory is closing and we are back to Yukawa's meson theory, except that we 
have learned to add one important refinement to the theory: broken chiral 
symmetry is a crucial constraint that generates and controls the dynamics 
and establishes a clear connection with the underlying theory, QCD. 

Following the first initiative by Weinberg [18], pioneering work was per- 
formed by Ordonez, Ray, and van Kolck |1 9|. 120] who constructed a 7V./V 
potential in coordinate space based upon chiral perturbation theory at next- 
to-next-to-leading order. The results were encouraging and many researchers 
became attracted to the new field [SU EH [23l EH [251 ESI E7] . As a conse- 
quence, nuclear EFT has developed into one of the most popular branches 
of modern nuclear physics |28[ [29] . 

It is the purpose of these lectures to describe in some detail the recent 
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progress in our understanding of nuclear forces in terms of nuclear EFT. 

2 QCD and the Nuclear Force 

Quantum chromo dynamics (QCD) is the theory of strong interactions. It 
deals with quarks, gluons and their interactions and is part of the Standard 
Model of Particle Physics. QCD is a non-Abelian gauge field theory with 
color SU(3) the underlying gauge group. The non-Abelian nature of the 
theory has dramatic consequences. While the interaction between colored 
objects is weak at short distances or high momentum transfer ( "asymptotic 
freedom"); it is strong at long distances ( ~ 1 fm) or low energies, leading 
to the confinement of quarks into colorless objects, the hadrons. Conse- 
quently, QCD allows for a perturbative analysis at large energies, whereas it 
is highly non-perturbative in the low-energy regime. Nuclear physics resides 
at low energies and the force between nucleons is a residual QCD interac- 
tion. Therefore, in terms of quarks and gluons, the nuclear force is a very 
complicated problem. 

3 Effective Field Theory for Low-Energy QCD 

The way out of the dilemma of how to derive the nuclear force from QCD 
is provided by the effective field theory (EFT) concept. First, one needs to 
identify the relevant degrees of freedom. For the ground state and the low- 
energy excitation spectrum of an atomic nucleus as well as for conventional 
nuclear reactions, quarks and gluons are ineffective degrees of freedom, while 
nucleons and pions are the appropriate ones. Second; to make sure that this 
EFT is not just another phenomenology, the EFT must observe all relevant 
symmetries of the underlying theory. This requirement is based upon a 'folk 
theorem' by Weinberg |17j : 

If one writes down the most general possible Lagrangian, in- 
cluding all terms consistent with assumed symmetry principles, 
and then calculates matrix elements with this Lagrangian to any 
given order of perturbation theory, the result will simply be the 
most general possible S-matrix consistent with analyticity, per- 
turbative unitarity, cluster decomposition, and the assumed sym- 
metry principles. 

Thus, the EFT program consists of the following steps: 



5 



1 . Identify the degrees of freedom relevant at the resolution scale of (low- 
energy) nuclear physics: nucleons and pions. 

2. Identify the relevant symmetries of low-energy QCD and investigate if 
and how they are broken. 

3. Construct the most general Lagrangian consistent with those symme- 
tries and the symmetry breaking. 

4. Design an organizational scheme that can distinguish between more 
and less important contributions: a low-momentum expansion. 

5. Guided by the expansion, calculate Feynman diagrams to the the de- 
sired accuracy for the problem under consideration. 

We will now elaborate on these steps, one by one. 
3.1 Symmetries of Low- Energy QCD 

In this section, we will give a brief introduction into (low-energy) QCD, its 
symmetries and symmetry breaking. A more detailed introduction can be 
found in the excellent lecture series by Scherer and Schindler [30 . 

3.1.1 Chiral Symmetry 

The QCD Lagrangian reads 

£ Q CD = q(il^ - M)q - \Q^ V (1) 
with the gauge-covariant derivative 

T>p = d„ + ig^-A^a (2) 
and the gluon field strength tensor 

In the above, q denotes the quark fields and M. the quark mass matrix. 
Further, g is the strong coupling constant and A^^ a are the gluon fields. 
The A a are the Gell-Mann matrices and the f a b c the structure constants 
of the S , C/(3) co ior Lie algebra (a,b,c = 1, . . . , 8); summation over repeated 
indices is always implied. The gluon-gluon term in the last equation arises 
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from the non-Abelian nature of the gauge theory and is the reason for the 
peculiar features of the color force. 

On a typical hadronic scale, i.e., on a scale of low-mass hadrons which are 
not Goldstone bosons, e.g., m p = 0.78 GeV « 1 GeV; the masses of the up 
(u), down (d), and — to a certain extend — strange (s) quarks are small 



m u = 2 ± 1 MeV (4) 
m d = 5 ± 2 MeV (5) 
m s = 95 ± 25 MeV (6) 

It is therefore of interest to discuss the QCD Lagrangian in the limit of 
vanishing quark masses: 

4cD = ^7^-^WSr- (7) 
Defining right- and left-handed quark fields, 

QR = PrQ , 1L = Pm , (8) 

with 

^ = i(l + 75), ^ = ^(1-75), (9) 
we can rewrite the Lagrangian as follows: 

£qcd = QRH^qR + QLvfVpqL ~ \G^aG» v ■ (10) 

Restricting ourselves now to up and down quarks, we see that £qcd i s 
invariant under the global unitary transformations 



and 



[I )[l • ^ 



where n (i = 1,2,3) are the generators of 5f7(2)fl avor , the usual Pauli spin 
matrices. The right- and left-handed components of massless quarks do not 
mix. This is SU(2)r x SU(2)r symmetry, also known as chiral symmetry. 
Noether's Theorem implies the existence of six conserved currents; three 
right-handed currents 

Ri = QR1*~<1R with <9 M i?f = (13) 
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and three left-handed currents 

Lt = qLl^q L with 8^ = 0. (14) 

It is useful to consider the following linear combinations; namely, three vec- 
tor currents 

Vf = R? + L? = q^q with 8^ = (15) 

and three axial-vector currents 

A? = R?-L? = 9 7 "75^Z with 8^ = 0, (16) 

which got their names from the fact that they transform as vectors and 
axial- vectors, respectively Thus, the chiral SU(2)i x SU(2)r symmetry is 
equivalent to SU{2)y x SU(2)a, where the vector and axial-vector transfor- 
mations are given respectively by 

and 

q = ( d J ' ~^ V ,0 2 J \ d f 

Obviously, the vector transformations are isospin rotations and, therefore, 
invariance under vector transformations can be identified with isospin sym- 
metry. 

There are the six conserved charges, 

v 



' -exp -^75? " . (18) 



QY = J d 3 x V? = J d 3 x q\t,x)^q{t,x) with = (19) 



dQ 

= / a x v{ = j a~x q~ {i,x) — q{i,x) witn 
and 

Qt = Jd 3 xA° = J d 3 x q\t, x) l5 ^q(t, x) with = , (20) 

which are also generators of SU(2)y x SU(2)a- 
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3.1.2 Explicit Symmetry Breaking 

The mass term —qAiq in the QCD Lagrangian Eq. ([I]) breaks chiral sym- 
metry explicitly. To better see this, let's rewrite Ai, 



m = 7 : (2D 





(m u + m d )[ * " )+^{m u -m d )[ " " I (22) 



2 

= ^(m u + m d ) I + -(m u - m d ) t 3 . (23) 

The first term in the last equation in invariant under SU(2)y (isospin sym- 
metry) and the second term vanishes for m u = m d . Thus, isospin is an exact 



symmetry if m u = m d . However, both terms in Eq. (23) break SU(2)a- 
Since the up and down quark masses are small as compared to the typical 
hadronic mass scale of ~ 1 GeV [cf. Eqs. Q and (|5j], the explicit chiral 
symmetry breaking due to non- vanishing quark masses is very small. 

3.1.3 Spontaneous Symmetry Breaking 

A (continuous) symmetry is said to be spontaneously broken if a symmetry 
of the Lagrangian is not realized in the ground state of the system. There is 
evidence that the chiral symmetry of the QCD Lagrangian is spontaneously 
broken — for dynamical reasons of nonperturbative origin which are not fully 
understood at this time. The most plausible evidence comes from the hadron 
spectrum. From chiral symmetry, one would naively expect the existence 
of degenerate hadron multiplets of opposite parity, i.e., for any hadron of 
positive parity one would expect a degenerate hadron state of negative parity 
and vice versa. However, these "parity doublets" are not observed in nature. 
For example, take the />meson, a vector meson with negative parity (1 — ) 
and mass 776 MeV. There does exist a 1 + meson, the ai, but it has a mass 
of 1230 MeV and, thus, cannot be perceived as degenerate with the p. On 
the other hand, the p meson comes in three charge states (equivalent to 
three isospin states), the p^- and the p° with masses that differ by at most 
a few MeV. In summary, in the QCD ground state (the hadron spectrum) 
SU(2)v (isospin symmetry) is well observed, while SU(2)a (axial symmetry) 
is broken. Or, in other words, SU(2)y x SU(2)a is broken down to SU(2)v- 
A spontaneously broken global symmetry implies the existence of (mass- 
less) Goldstone bosons with the quantum numbers of the broken generators. 



The broken generators are the Qf of Eq. (20) which are pseudoscalar. The 
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Goldstone bosons are identified with the isospin triplet of the (pseudoscalar) 
pions, which explains why pions are so light. The pion masses are not ex- 
actly zero because the up and down quark masses are not exactly zero either 
(explicit symmetry breaking). Thus, pions are a truly remarkable species: 
they reflect spontaneous as well as explicit symmetry breaking. 



3.2 Chiral Effective Lagrangians Involving Pions 

The next step in our EFT program is to build the most general Lagrangian 
consistent with the (broken) symmetries discussed above. An elegant formal- 
ism for the construction of such Lagrangians was developed by Callan, Cole- 
man, Wess, and Zumino (CCWZ) [32 who worked out the group-theoretical 
foundations of non-linear realizations of chiral symmetry. The Lagrangians 
given below are built upon the CCWZ formalism. 

As discussed, the relevant degrees of freedom are pions (Goldstone bosons) 
and nucleons. Since the interactions of Goldstone bosons must vanish at zero 
momentum transfer and in the chiral limit (m — ► 0), the low-energy expan- 
sion of the Lagrangian is arranged in powers of derivatives and pion masses. 
This is chiral perturbation theory (ChPT). 

The Lagrangian consists of one part that deals with the interaction 
among pions, C nn , and another one that describes the interaction between 
pions and the nucleon, C n N- 

£eff = £tttv + £-kN (24) 

with 

£-inr = ^7T7r + ^-tttt + • ■ ■ (25) 

and 

£ 7rJV = 4 1 i + 4 2 i + 4 3 i + --- ) (26) 

where the superscript refers to the number of derivatives or pion mass inser- 
tions (chiral dimension) and the ellipsis stands for terms of higher dimension. 
The leading order (LO) irn Lagrangian is given by [33] 

C^ = f -fti[d^Ud^ + ml{U + U ] )] (27) 
and the LO relativistic ttN Lagrangian reads 

C<$ f = *(i'fD li -M N + y7 M ~ ->«,,) * (2*) 
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with 



D„ 



3 M + r M 



i 

1 



itfd ll S-Zd tl tf) = ~T-d li n + ... 



u 



r = 1 + -f-T • TV 

Jit 



1 



r7T 



ia , , 3 8a — 1 4 



f 3 
j - 



7T^ + 



(29) 
(30) 

(31) 

(32) 



In Eq. (28) the chirally covariant derivative Dn is applied which introduces 



the "g auge term" (also known as chiral connection), a vector current that 
leads to a coupling of pions with the nucleon. Besides this, the Lagrangian 
includes a coupling term which involves the axial vector Un- The SU(2) 
matrix U = £ 2 collects the Goldstone pion fields. 

In the above equations, Mn denotes the nucleon mass, qa the axial- 
vector coupling constant, and f n the pion decay constant. Numerical values 
will be given later. 



The coefficient a that appears in Eq. (32) is arbitrary. Therefore, dia- 



grams with chiral vertices that involve three or four pions must always be 
grouped together such that the a-dependence drops out (cf. Fig. [4] below). 

We apply the heavy baryon (HB) formulation of chiral perturbation the- 
ory [35: in which the relativistic ttN Lagrangian is subjected to an expansion 
in terms of powers of 1 /Mn (kind of a nonrelativistic expansion) , the lowest 
order of which is 



C 



(1) 

ttN 



N iD, 







— a ■ u iV 
2 



N 



id 



1 



7T 



x di 



'07TJ 



9A_ 
2U 



T ■ ((T ■ V)7T 



N + ... (33) 



In the relativistic formulation, the nucleon is represented by a four-component 
Dirac spinor field, ^f, while in the HB version, the nucleon, N, is a Pauli 
spinor; in addition, all nucleon fields include Pauli spinors describing the 
isospin of the nucleon. 

At dimension two, the relativistic ttN Lagrangian reads 



( 2 ), 



(34) 



8=1 



The various operators 0\ are given in Ref. [36]. The fundamental rule by 
which this Lagrangian — as well as all the other ones — are assembled is that 
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they must contain all terms consistent with chiral symmetry and Lorentz 
invariance (apart from other trivial symmetries) at a given chiral dimension 
(here: order two). The parameters Cj are known as low-energy constants 
(LECs) and are determined empirically from fits to ttN data. 

The HB projected irN Lagrangian at order two is most conveniently 
broken up into two pieces, 



with 



and 



£ 



?(2) 

-nN, ct 



(2) 

nN. fix 



N 



£ 



N 



(2) 
ttN 



c 



(2) 

7riV,fix 



+ £ (2) 

T 7rN, Ct ' 



2M N 



D-D + i 



9A 
4M N 



{a ■ D,u } 



N 



(35) 
(36) 



i 

+ 2 



2 ci ml (U + tf) + ^c 2 - 

( C4 + ik) d " ( " x " ) 



8Mtv y 



^0 + C 3 U H U 



N. 



(37) 



Note that C L fix is created entirely from the HB expansion of the relativistic 



jC^y and thus has no free parameters ( "fixed" ) , while £^y ct is dominated 
by the new ttN contact terms proportional to the Cj parameters, besides 
some small 1/Mjy corrections. 

At dimension three, the relativistic irN Lagrangian can be formally writ- 
ten as 



(2) 



2:-! 



(3), 



(38) 



with the operators, o\ 3 \ listed in Refs. (SHI EZ]; not all 23 terms are of 
interest here. The new LECs that occur at this order are the di. Similar 
to the order two case, the HB projected Lagrangian at order three can be 
broken into two pieces, 



£ 



(3) 

ttN 



(3) 

7rAf,fix 



+ £ (3) 

T '-'nN, ct ' 



(39) 



with £^ fix and £^ ct given in Refs. [361 EH- 



5(3) 



3.3 Nucleon Contact Lagrangians 

Nucleon contact interactions consist of four nucleon fields (four nucleon legs) 
and no meson fields. Such terms are needed to renormalize loop integrals, 
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to make results reasonably independent of regulators, and to parametrize 
the unresolved short-distance contributions to the nuclear force. For more 



about contact terms, see Sec. 5.3 



Because of parity, nucleon contact interactions come only in even num- 
bers of derivatives, thus, 

CNN = + £ { NN + £ { N N + ■ ■ ■ (40) 

The lowest order (or leading order) NN Lagrangian has no derivatives 
and reads [T8] 

C$ N = -^C S NNNN - -C T NaNNaN , (41) 

where N is the heavy baryon nucleon field. Cs and Ct are unknown con- 
stants which are determined by a fit to the NN data. The second order NN 
Lagrangian is given by |19] 

C$ N = -C[[(NVN) 2 + (VNN) 2 ]-C' 2 {NVN)-(VNN) 
-C' 3 NN[NV 2 N + V 2 iViV] 

-iC'^NVN • (VJV x ffN) + (VN)N ■ (Na x VN)\ 
-iC' 5 NN(VN ■ a x VN) - iC' 6 (NaN) ■ (VN x ViV) 
— (C'jdikSji + C' 8 5u5kj + C' 9 5ijSki) 
x [Na k diNNaidjN + d^Na k NdjNaiN] 



-{C'^ikSji + C n Sa5 k j + C'^SijSkONakdiNdjNaiN 



+C u 8 ij 6ki) [diNakdjN + djNakdi^NaiN . (42) 

Similar to Cs and Ct, the C[ are unknown constants which are fixed in a 
fit to the iViV data. Obviously, these contact Lagrangians blow up quite a 
bit with increasing order, which why we do not give C^N explicitly here. 



4 Nuclear Forces from EFT: Overview 

In the beginning of Sec. [3j we spelled out the steps we have to take to 
accomplish our EFT program for the derivation of nuclear forces. So far, 
we discussed steps one to three. What is left are steps four (low-momentum 
expansion) and five (Feynman diagrams). In this section, we will say more 
about the expansion we are using and give an overview of the Feynman 
diagrams that arise order by order. 
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4.1 Chiral Perturbation Theory and Power Counting 

In ChPT, we analyze contributions in terms of powers of small momenta 
over the large scale: (Q/A x ) u , where Q stands for a momentum (nucleon 
three-momentum or pion four-momentum) or a pion mass and A x ~ 1 GeV 
is the chiral symmetry breaking scale (hadronic scale). Determining the 
power v at which a given diagram contributes has become known as power 
counting. For a non-iterative contribution involving A nucleons, the power 
v is given by 

v = -2 + 2 A - 2C + 2L + ^ A< , (43) 

i 

with 

+ (44) 

where C denotes the number of separately connected pieces and L the num- 
ber of loops in the diagram; di is the number of derivatives or pion-mass 
insertions and rtj the number of nucleon fields involved in vertex i; the sum 
runs over all vertices contained in the diagram under consideration. Note 
that for an irreducible NN diagram (A = 2), the above formula reduces to 

v = 2L + J2 A i ( 45 ) 

i 

The power v is bounded from below; e.g., for A = 2, v > 0. This fact is 
crucial for the power expansion to be of any use. 

4.2 The Hierarchy of Nuclear Forces 

Chiral perturbation theory and power counting imply that nuclear forces 
emerge as a hierarchy ruled by the power u, Fig. [T] 

The NN amplitude is determined by two classes of contributions: con- 
tact terms and pion-exchange diagrams. There are two contacts of order 
Q° [0(Q )] represented by the four-nucleon graph with a small-dot vertex 
shown in the first row of Fig. [TJ The corresponding graph in the second row, 
four nucleon legs and a solid square, represents the seven contact terms of 
0(Q 2 ). Finally, at 0(Q 4 ), we have 15 contact contributions represented by 
a four-nucleon graph with a solid diamond. 

Now, turning to the pion contributions: At leading order [LO, O(Q ), 
v = 0], there is only the well-known static one-pion exchange (1PE), second 
diagram in the first row of Fig. [T] Two-pion exchange (2PE) starts at next- 
to-leading order (NLO, v = 2) and all diagrams of this leading-order two- 
pion exchange are shown. Further 2PE contributions occur in any higher 
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2N Force 



3N Force 4N Force 



LO 



Cr 

NLO 



Q- 



NNLO 




Q 4 

N'LO 



Figure 1: Hierarchy of nuclear forces in ChPT. Solid lines represent nucleons 
and dashed lines pions. Further explanations are given in the text. 



order. Of this sub-leading 2PE, we show only two representative diagrams at 
next-to-next-to-leading order (NNLO) and three diagrams at next-to-next- 
to-next-to-leading order (N 3 LO). 

Finally, there is also three-pion exchange, which shows up for the first 
time at N 3 LO (two loops; one representative 3tt diagram is included in 
Fig. [T|. At this order, the Sir contribution is negligible [38]. 

One important advantage of ChPT is that it makes specific predictions 
also for many-body forces. For a given order of ChPT, two-nucleon forces 
(2NF), three-nucleon forces (3NF), . . . are generated on the same footing (cf. 
Fig. [TJ. At LO, there are no 3NF, and at NLO, all 3NF terms cancel (IS) 
[39] . However, at NNLO and higher orders, well-defined, nonvanishing 3NF 
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occur |39| 14*0] . Since 3NF show up for the first time at NNLO, they are 
weak. Four-nucleon forces (4NF) occur first at N 3 LO and, therefore, they 
are even weaker. 



5 Two-Nucleon Forces 

In this section, we will elaborate in detail on the two-nucleon force contri- 
butions of which we have given a rough overview in the previous section. 

5.1 Pion-Exchange Contributions in ChPT 



The effective pion Lagrangians presented in Sec. 3.2 are the crucial ingredi- 
ents for the evaluation of the pion-exchange contributions to the NN inter- 
action. We will derive these contributions now order by order. 

We will state our results in terms of contributions to the momentum- 
space NN amplitude in the center-of-mass system (CMS), which takes the 
general form 

V(p',p) = Vc +Tl-T 2 W C 

+ [V s + ri ■ r 2 W s } ai ■ a 2 

+ [V LS + T 1 -T 2 W L s\(-iS-(q>ik)) 

+ [Vr + Ti ■ t 2 Wt ] <J\- qo 2 - q 

+ [Vvl + ti ■ T 2 W aL ] ?i-(?x k) a 2 -(qxk), (46) 

where p ' and p denote the final and initial nucleon momenta in the CMS, 
respectively; moreover, 

q = p'—p is the momentum transfer, 

k = \{p ' + p) the average momentum, (47) 

S = + <?2) the total spin, 

and a\ t 2 and T\ j2 are the spin and isospin operators, respectively, of nucleon 
1 and 2. For on-energy-shell scattering, V a and W a (a = C, S, LS, T, aL) 
can be expressed as functions of q and k (with q = \q\ and k = \k\), only. 

Our formalism is similar to the one used by the Munich group |22| |4"T| |4"2] 
except for two differences: all our momentum space amplitudes differ by an 
over-all factor of (—1) and our spin-orbit potentials, Vls and Wis-, differ by 
an additional factor of (—2). Our conventions are more in tune with what 
is commonly used in nuclear physics. 
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In all expressions given below, we will state only the nonpolynomial con- 
tributions to the NN amplitude. Note, however, that dimensional regular- 
ization typically generates also polynomial terms. These polynomials are 
absorbed by the contact interactions to be discussed in a later section and, 
therefore, they are of no interest here. 



5.1.1 Zeroth Order (LO) 

At order zero [u = 0, O(Q ), lowest order, leading order, LO], there is only 
the well-known static one-pion exchange, second diagram in the first row of 
Fig. [T] which is given by: 



Vln(p ,P) = -T72 T l • T 2 —21 2~ ■ ( 48 

At first order [u = 1, 0(Q)], there are no pion-exchange contributions 
(and also no contact terms). 



5.1.2 Second Order (NLO) 

Non- vanishing higher-order graphs start at second order {y = 2, next-to- 
leading order, NLO). The most efficient way to evaluate these loop diagrams 
is to use covariant perturbation theory and dimensional regularization. This 
is the method applied by the Munich group |22| W\\ 142] . One starts with 



the relativistic versions of the irN Lagrangians (cf. Sec. 3.2 ) and sets up 
four-dimensional (covariant) loop integrals. Relativistic vertices and nucleon 
propagators are then expanded in powers of 1/M/v- The divergences that 
occur in conjunction with the four-dimensional loop integrals are treated by 
means of dimensional regularization, a prescription which is consistent with 
chiral symmetry and power counting. The results derived in this way are 
the same obtained when starting right away with the HB versions of the 
ttN Lagrangians. However, as it turns out, the method used by the Munich 
group is more efficient in dealing with the rather tedious calculations. 

Two-pion exchange occurs first at second order, also know as leading- 
order 2tt exchange. The graphs are shown in the first row of Fig. [2] Since 
a loop creates already v = 2, the vertices involved at this order can only be 
from the leading/lowest order Lagrangian C^jf, Eq. (33), i. e., they carry 



only one derivative. These vertices are denoted by small dots in Fig. [2] 
Concerning the box diagram, we should note that we include only the non- 
iterative part of this diagram which is obtained by subtracting the iter- 



ated 1PE contribution Eq. (65) or Eq. (66), below, but using M^/E. 
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Q 2 

(NLO) 



Q 3 

(r^LO) 



Figure 2: Twopion exchange contributions to the NN interaction at order two 
and three in small momenta. Solid lines represent nucleons and dashed lines p ions . 



Small dots denote vertices from the leading order irN Lagrangian Eq. (133 > 



Large solid dots are vertices proportional to the LECs c, from the second order 



Lagrangian A^jy ct , Eq. (37 1. Symbols with an open circles are relativistic 1/Mjv 
corrections contained in the second order Lagrangian 

Eqs. 135b. Only a few 
representative examples of 1/Mjv corrections are shown and not all. 



M 2 N /E p n M N at this order (NLO). Summarizing all contributions from 
irreducible two-pion exchange at second order, one obtains |22| : 



W c 



V T 



where 



and 



+ 



m 



4ml(5g A A - Ag\ - 1) + q 2 {22,g\ - 10<& - 1) 



647T2/4 ' 

m 



w . w + q 

— In 

q 2m^ 



w 



Am 2 + q 2 . 



(49) 
(50) 

(51) 
(52) 
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5.1.3 Third Order (NNLO) 

The two-pion exchange diagrams of order three (y = 3, next-to-next-to- 
leading order, NNLO) are very similar to the ones of order two, except that 
they contain one insertion from Eq. (351). The resulting contributions 

are typically either proportional to one of the low-energy constants Cj or 
they contain a factor 1/Mjy. Notice that relativistic 1/Mjv corrections can 
occur for vertices and nucleon propagators. In Fig. [2j we show in row 2 
the diagrams with vertices proportional to q (large solid dot), Eq. (37), 
and in row 3 and 4 a few representative graphs with a l/M^r correction 
(symbols with an open circle). The number of l/M^v correction graphs 
is large and not all are shown in the figure. Again, the box diagram is 
corrected for a contribution from the iterated IPE. If the iterative 2PE of 
Eq. (65) is used, the expansion of the factor M^/E p = Mn — p 2 /2Mn + • • • 



is applied and the term proportional to (—p 2 /2M]y) is subtracted from the 
third order box diagram contribution. Then, one obtains for the full third 
order contribution |22|: 



V c 
W c 

V T 
W t 



39a j 9 A m l 



16vr/4 1 16M N w 2 



2m 2 n (2ci - c 3 ) -q 2 \o i + 



16M N J 



x w 



91 



{39A m l w ~ 



\2MM N f± 
- 'Ami + ~ 9A^rnl + 3q 2 ) 
9g\w 2 A(q) 



w 



l A(q)} 



V LS 
W ls 
with 



-\w s 

gZ 

9\A{q) 
327T/4 

39\w 2 A{q) 

gj(! -9a) . .2 

327TMjv/4 



5127rMjv/4 



C4 + 



4M N 



w 



_A_ 
8M N 



(10ml + 3g 2 ) 



■ 4 w 2 A(q) , 



1 



A(q) = — arctan 

2q 2m„ 



(53) 

(54) 
(55) 

(56) 
(57) 
(58) 

(59) 



19 



and 



w 



2m\ + I 2 • 



(60) 



As discussed in Sec. 5.1.5 below, we prefer the iterative 2PE defined in 
(66), which leads to a different NNLO term for the iterative 2PE. This 
changes the 1/Mjy terms in the above potentials. The changes are obtained 



Eq 



by adding to Eqs. (53)-(56) the following terms: 



v - ~ -7- Vs = ^M^ im ' +w2AW,) <63) 

WT - -7 Ws '-vJfiM-« (m * +u2m) <64) 



5.1.4 Fourth Order (N 3 LO) 

This order, which may also be denoted by next-to-next-to-next-to- leading 
order (N 3 LO), is very involved. Three-pion exchange (3PE) occurs for the 
first time at this order. The 3PE contribution at N 3 LO has been calculated 
by the Munich group and found to be negligible [38]. Therefore, we will 
ignore it. 

The 2PE contributions at N 3 LO can be subdivided into two groups, one- 
loop graphs, Fig. [3] and two-loop diagrams, Fig. [4]. Since these contributions 
are very complicated, we have moved them to Appendix A. 



5.1.5 Iterated One-Pion- Exchange 

Besides all the irreducible 2PE contributions presented above, there is also 
the reducible 2PE which is generated from iterated 1PE. This "iterative 
2PE" is the only 2PE contribution which produces an imaginary part. Thus, 
one wishes to formulate this contribution such that relativistic elastic uni- 
tarity is satisfied. There are several ways to achieve this. 

Kaiser et al. [22 define the iterative 2PE contribution as follows, 



tt(kbw) t-> i -,_M% f dV 
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of 



(N 3 LO) 




Figure 3: One- loop 27r-exchange contributions to the NN interaction at order 
four. Basic notation as in Fig. [2j Symbols with a large solid dot and an open circle 
denote 1/Mjv corrections of vertices proportional to Cj. Symbols with two open 
circles mark relativistic 1/ Mfj corrections. Both corrections are part of the third 
order Lagrangian Eq. ( |39[ ). Representative examples for all types of one- loop 

graphs that occur at this order are shown. 



with Vi n given in Eq. (48) 



Since we adopt the relativistic scheme developed by Blankenbecler and 
Sugar jl3] (BbS) (see beginning of Sec. 5.4), we prefer the following for- 



mulation which is consistent with the BbS approach (and, of course, with 
relativistic elastic unitarity): 



Vr 



(EM) 



2nAt 



{p',p) 



d 3 p" M 2 N V ln (p',p")V l7T (p",p) 



(2vr) 3 E p „ 



P 



p" 2 + ie 



(66) 



The iterative 2PE contribution has to be subtracted from the covariant 
box diagram, order by order. For this, the expansion M^/E p = Mn — 
p 2 /2M N + ... is applied in Eq. ^ and M 2 N jE v „ = M N -p" 2 /2M N + ... in 
Eq. ((661). At NLO, both choices for the iterative 2PE collapse to the same, 
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Q 4 

(N 3 LO) 




Figure 4: Two-loop 27r-exchange contributions at order four. Basic notation as in 
Fig. [2] The oval stands for all one-loop irN graphs some of which are shown in the 
lower part of the figure. The solid square represents vertices proporti onal to the 
LECs di which are introduced by the third order Lagrangian £^y, Eq. (38 1. More 
explanations are given in the text. 



while at NNLO there are obvious differences. 

5.2 NN Scattering in Peripheral Partial Waves Using the 
Perturbative Amplitude 

After the tedious mathematics of the previous section, it is time for more 
tangible affairs. The obvious question to address now is: How does the 
derived NN amplitude compare to empirical information? Since our deriva- 
tion includes only one- and two-pion exchanges, we are dealing here with 
the long- and intermediate-range part of the NN interaction. This part of 
the nuclear force is probed in the peripheral partial waves of NN scattering. 
Thus, in this section, we will calculate the phase shifts that result from the 
NN amplitudes presented in the previous section and compare them to the 
empirical phase shifts as well as to the predictions from conventional meson 
theory. Besides the irreducible two-pion exchanges derived above, we must 
also include 1PE and iterated 1PE. 

In this section |44j . which is restricted to just peripheral waves, we 
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will always consider neutron-proton (np) scattering and take the charge- 
dependence of 1PE due to pion-mass splitting into account, since it is ap- 
preciable. With the definition 

4f4 q 2 + mi 
the charge-dependent 1PE for np scattering is 

V{7\p',p) = + ("l) m 214K ± ) , (68) 

where I denotes the isospin of the two-nucleon system. We use m w o = 
134.9766 MeV, m n ± = 139.5702 MeV [M], and 

M N = ™ pM Jl = 938.9182 MeV . (69) 
M p + M n v ' 

Also in the iterative 2PE, we apply the charge-dependent 1PE, i.e., in 



Eq. ( 66 ) we replace V\ v with V"/ . 



The perturbative relativistic T-matrix for np scattering in peripheral 
waves is 

T(p',p) = V}: p \p',p) + V^ np \p',p) + V^, irr (p',P) » (70) 
where V^-Kirr refers to any or all of the irreducible 2PE contributions pre- 



sented in Sec. 5.1 depending on the order at which the calculation is con- 
ducted. In the calculation of the irreducible 2PE, we use the average pion 
mass m n = 138.039 MeV and, thus, neglect the charge-dependence due to 
pion-mass splitting. The charge-dependence that emerges from irreducible 
27T exchange was investigated in Ref. jl5] and found to be negligible for 
partial waves with L > 3. 



For the T-matrix given in Eq. ( 70 ) , we calculate phase shifts for partial 



waves with L > 3 and Tiaj, < 300 MeV. At order four in small momenta, 
partial waves with L > 3 do not receive any contributions from contact inter- 
actions and, thus, the non-polynomial pion contributions uniquely predict 
the F and higher partial waves. We use f n = 92.4 MeV [31] and gA = 1.29. 
Via the Goldberger-Treiman relation, gA = g-wNN fn/M^, our value for gA 
is consistent with g^Niy/^ir = 13.63 ± 0.20 which is obtained from 7riV and 
NN analysis HUE]. 

The LECs used in this calculation are shown in Table [2] column "NN 
periph. Fig. Note that many determinations of the LECs, Cj and di, can 
be found in the literature. The most reliable way to determine the LECs 
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Table 2: Low-energy constants, LECs, used for a NN potential at N 3 LO, 
Sec. |5.4.4| and in the calculation of the peripheral NN phase shifts shown 
in Fig. [5] (column "NN periph. Fig. pjj'). The Cj belong to the dimension- 
two 7tN Lagrangian, Eq. (37), and are in units of GeV -1 , while the di are 



associated with the dimension-three Lagrangian, Eq. (38 ), and are in units of 
~ 2 The column u ttN empirical" shows determinations from irN data. 



GeV" 



LEC 


NN potential 


NN periph. 


vriV 




at N 3 LO 


Fig.g 


empirical 


Cl 


-0.81 


-0.81 


-0.81 ±0.15 a 


C2 


2.80 


3.28 


3.28 ± 0.23 6 


C3 


-3.20 


-3.40 


-4.69 ± 1.34 a 


c 4 


5.40 


3.40 


3.40 ± 0.04 a 


d\ + d 2 


3.06 


3.06 


3.06±0.21 6 


ds 


-3.27 


-3.27 


-3.27±0.73 fe 


d 5 


0.45 


0.45 


0.45 ± 0.42 6 


du — di*, 


-5.65 


-5.65 


-5.65 ±0.41 6 


"Table 1, 


Fit 1 of Ref. g8]. 







b Table 2, Fit 1 of Ref. [37]. 



from empirical irN information is to extract them from the ttN amplitude 
inside the Mandelstam triangle (unphysical region) which can be constructed 
with the help of dispersion relations from empirical irN data. This method 
was used by Biittiker and Meifiner [48 j . Unfortunately, the values for C2 and 
all di parameters obtained in Ref. [48 a carry uncertainties, so large that the 
values cannot provide any guidance. Therefore, in Table [2j only ci, C3, and 
C4 are from Ref. jH], while the other LECs are taken from Ref. [37] where 
the 7riV amplitude in the physical region was considered. To establish a link 
between ttN and iV./V, we apply the values from the above determinations 
in our calculations of the NN peripheral phase shifts. In general, we use 
the mean values; the only exception is C3, where we choose a value that is, 
in terms of magnitude, about one standard deviation below the one from 
Ref. [H]. With the exception of C3, phase shift predictions do not depend 
sensitively on variations of the LECs within the quoted uncertainties. 

In Fig. [5] we show the phase-shift predictions for neutron-proton scat- 
tering in F waves for laboratory kinetic energies below 300 MeV (for G and 
H waves, see Ref. [26 ). The orders displayed are defined as follows: 



• Leading order (LO) is just 1PE, Eq. (68). 
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Figure 5: F-wave phase shifts of neutron-proton scattering for laboratory kinetic 
energies below 300 MeV. We show the predictions from chiral pion exchange at lead- 
ing order (LO), next-to-leading order (NLO), next-to-next-to-leading order (N2LO), 
and next-to-next-to-next-to-leading order (N3LO). The solid dots and open circles 
are the results from the Nijmegen multi-energy np phase shift analysis |49j and the 
VPI single-energy np analysis SM99 |SD], respectively. 



Next-to-leading order (NLO) is 1PE, Eq. ([68]), plus iterated 1PE 



Eq. (66), plus the contributions of Sec. 5.1.2 (order two), Eqs. (49) 
and (50). 



Next-to-next-to-leading order (denoted by N2LO in the figures) con- 
sists of NLO plus the contributions of Sec. 5.1.3| (order three), Eqs. (53) 
d58b and (|6lb-(|64b. 



Next-to-next-to-next-to-leading order (denoted by N3LO in the fig- 



ures) consists of N2LO plus the contributions of Sec. 5.1.4 (order four), 
Eqs. (f99])-(pJL2|) and (flT5]-(p4l). 



It is clearly seen in Fig. [5] that the leading order 2-7T exchange (NLO) is 
a rather small contribution, insufficient to explain the empirical facts. In 
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Figure 6: F-wave phase shifts of neutron-proton scattering for laboratory kinetic 
energies below 300 MeV. We show the results from one-pion-exchange (OPE), 
and one- plus two-pion exchange as predicted by ChPT at next-to-next-to-next- 
to-leading order (N3LO) and by the Bonn Full Model [13] (Bonn). Note that the 
"Bonn" curve does not include the repulsive to and 7rp exchanges of the full model, 
since this figure serves the purpose to compare just predictions by different mod- 
els/theories for the 7r + 2ir contribution to the NN interaction. Empirical phase 
shifts (solid dots and open circles) as in Fig. [5] 



contrast, the next order (N2LO) is very large, several times NLO. This is 
due to the irirNN contact interactions proportional to the LECs Cj that are 
introduced by the second order Lagrangian C n jy, Eq. (34). These contacts 
are supposed to simulate the contributions from intermediate A-isobars and 
correlated 2tt exchange which are known to be large (see, e. g., Ref. |13j). 

At N3LO a clearly identifiable trend towards convergence emerges. Ob- 
viously, F$ and 3 F4 appear fully converged. However, in 3 i<2 and 3 i<3, N3LO 
differs noticeably from NNLO, but the difference is much smaller than the 
one between NNLO and NLO. This is what we perceive as a trend towards 
convergence. 
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In Fig. |6j we conduct a comparison between the predictions from chi- 
ral one- and two-pion exchange at N3LO and the corresponding predictions 
from conventional meson theory (curve 'Bonn'). As representative for con- 
ventional meson theory, we choose the Bonn meson-exchange model for the 
NN interaction [13], since it contains a comprehensive and thoughtfully con- 
structed model for 2n exchange. This 2ir model includes box and crossed 
box diagrams with NN, NA, and AA intermediate states as well as di- 
rect tttt interaction in S- and P-waves (of the tttt system) consistent with 
empirical information from irN and tttt scattering. Besides this the Bonn 
model also includes (repulsive) w-meson exchange and irreducible diagrams 
of 7r and p exchange (which are also repulsive). However, note that in the 
phase shift predictions displayed in Fig. [6j the "Bonn" curve includes only 
the i7r and 2tt contributions from the Bonn model; the short-range contri- 
butions are left out since the purpose of the figure is to compare different 
models/theories for it + 2n. In all waves shown we see, in general, good 
agreement between N3LO and Bonn. In 3 i*2 and 3 Fs above 150 MeV and 
in 3 i<4 above 250 MeV the chiral model at N3LO is more attractive than 
the Bonn 2tt model. Note, however, that the Bonn model is relativistic and, 
thus, includes relativistic corrections up to infinite orders. Thus, one may 
speculate that higher orders in ChPT may create some repulsion, moving 
the Bonn and the chiral predictions even closer together |51| . 

The 2tt exchange contribution to the NN interaction can also be de- 
rived from empirical ttN and tttt input using dispersion theory, which is 
based upon unitarity, causality (analyticity), and crossing symmetry. The 
amplitude NN — > tttt is constructed from ttN — > ttN and irN — > ttttN 
data using crossing properties and analytic continuation; this amplitude is 
then 'squared' to yield the NN amplitude which is related to NN by cross- 
ing symmetry [52]. The Paris group |11[ [T2] pursued this path and calcu- 
lated NN phase shifts in peripheral partial waves. Naively, the dispersion- 
theoretic approach is the ideal one, since it is based exclusively on empirical 
information. Unfortunately, in practice, quite a few uncertainties enter into 
the approach. First, there are ambiguities in the analytic continuation and, 
second, the dispersion integrals have to be cut off at a certain momentum 
to ensure reasonable results. In Ref. [13 a , a thorough comparison was con- 
ducted between the predictions by the Bonn model and the Paris approach 
and it was demonstrated that the Bonn predictions always lie comfortably 
within the range of uncertainty of the dispersion-theoretic results. There- 
fore, there is no need to perform a separate comparison of our chiral N3LO 
predictions with dispersion theory, since it would not add anything that we 
cannot conclude from Fig. [6] 
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Finally, we need to compare the predictions with the empirical phase 
shifts. In F waves the N 3 LO predictions above 200 MeV are, in general, 
too attractive. Note, however, that this is also true for the predictions by 
the Bonn tt + 2tt model. In the full Bonn model, besides tt + 2tt, (repul- 
sive) oj and up exchanges are included which move the predictions right 
on top of the data. The exchange of a u meson or combined irp exchange 
are 3tt exchanges. Three-pion exchange occurs first at chiral order four. 
It has be investigated by Kaiser [38] and found to be negligible, at this 
order. However, 3tt exchange at order five appears to be sizable [53] and 
may have impact on F waves. Besides this, there is the usual short-range 
phenomenology. In ChPT, this short-range interaction is parametrized in 
terms of four-nucleon contact terms (since heavy mesons do not have a place 
in that theory). Contact terms of order four (N 3 LO) do not contribute to 
_F-waves, but order six does. In summary, the remaining small discrepan- 
cies between the N3LO predictions and the empirical phase shifts may be 
straightened out in fifth or sixth order of ChPT. 



5.3 NN Contact Potentials 

In conventional meson theory, the short-range nuclear force is described by 
the exchange of heavy mesons, notably the u(782). The qualitative short- 
distance behavior of the iViV potential is obtained by Fourier transform of 
the propagator of a heavy meson, 

d 3 q^-^ ~ . (71) 

ChPT is an expansion in small momenta Q, too small to resolve struc- 
tures like a p(770) or w(782) meson, because Q <C A x m PtU . But the 
latter relation allows us to expand the propagator of a heavy meson into a 
power series, 

1 if Q 2 Q A \ 



m t + Q 2 m l)\ m t m t j 

where the uj is representative for any heavy meson of interest. The above 
expansion suggests that it should be possible to describe the short distance 
part of the nuclear force simply in terms of powers of Q/m w , which fits in 
well with our over- all power scheme since Q/m u ~ Q/A x . 

A second purpose of contact terms is renormalization. Dimensional reg- 



ularization of the loop integrals of pion-exchanges (cf. Sec. 5.1) typically 



2S 



generates polynomial terms with coefficients that are, in part, infinite or 
scale dependent. Contact terms pick up infinities and remove scale depen- 
dence. 

The partial-wave decomposition of a power Q v has an interesting prop- 
erty. First note that Q can only be either the momentum transfer between 



the two interacting nucleons q or the average momentum k [cf . Eq. ( 47 ) for 
their definitions]. In any case, for even u, 

Q» = f»(cos9), (73) 

where f m stands for a polynomial of degree m and 9 is the CMS scattering 
angle. The partial-wave decomposition of Q v for a state of orbital-angular 
momentum L involves the integral 



M 

1 L 



r+1 f+1 

J Q u P L (cos9)dcos9 = J f^(cos9)P L (cos9)dcos9, (74) 

where Pl is a Legendre polynomial. Due to the orthogonality of the Pl, 

I L v) = for L > V - . (75) 

Consequently, contact terms of order zero contribute only in S-waves, while 
order-two terms contribute up to P-waves, order-four terms up to D-waves, 
etc.. 

We will now present, one by one, the various orders of NN contact 
terms together with their partial- wave decomposition [SJj. Note that, due 
to parity, only even powers of Q are allowed. 

5.3.1 Zeroth Order 

The contact potential at order zero reads: 

V^(p',p) = C s + C T a 1 -a 2 (76) 

Partial wave decomposition yields: 

V^CS ) = a So =47r(C s -3C T ) 

V^{ 3 Si) = Cz Si =Ak(C s + C t ) (77) 
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5.3.2 Second Order 

The contact potential contribution of order two is given by: 

V^(p',p) = C iq 2 + C 2 k 2 

+ (c 3 q 2 + C A k 2 ) <7i • B 2 

+ C b (-iS-{qxk)) 
+ C Q {ai -q)(a 2 -q) 

+ C 7 (<ri ■ k) (a 2 ■ k) (78) 
Second order partial wave contributions: 

V^CSo) = C 1So (p 2 +p' 2 ) 

= 4vr (Ci + -Jc 2 - 3C7 3 - - C 6 - ^C 7 ) (p 2 + p' 2 ) 
1/ (2) ( 3 P ) = Cs PoP p' 

= 4rr + ^C 2 - |c 3 + \c, - + 2C 6 - lev) PP 

^CPi) = Ci Pl pp' 

= 4rr + ^C 2 + 2C 3 - + - ^C 7 ) pp' 

^ (2) ( 3 Si) = C 3Sl (p 2 +p' 2 ) 

= 4rr (c t + X -C2 + C 3 + \d + ^C 6 + ^C 7 ) (p 2 + p' 2 ) 

V (a) ( s 5i- 3 £>i) = C 3Sl _ 3Dl p 2 

= 4rr (- — C 6 -— C 7 jp 

V m ( 3 P 2 ) = C 3p2 pp' 

= 4tt (- |ci + ^C 2 - |c 3 + ^C 4 + ^C 5 ) PP (79) 

5.3.3 Fourth Order 

The contact potential contribution of order four reads: 

V^\p',p) = D iq 4 + D 2 k 4 + D 3 q 2 k 2 + D 4 (qxk) 2 

+ (D b q A + D 6 k 4 + D 7 q 2 k 2 + L> 8 (<? x A?) 2 ) «?! • a 2 
+ (iV + ^Diofc 2 ) • {q x fc)) 
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+ (D xl q 2 + D l2 k 2 ) (<?! • q) (a 2 ■ 9) 
+ (D m 2 + D u k 2 ) {B x ■ k) (a 2 ■ k) 
+ L>i5 (Vi • (9 X £) a 2 -(qxkj) (80) 

The rather lengthy partial- wave expressions of this order have been relegated 
to Appendix B. 

5.4 Constructing a Chiral iViV Potential 
5.4.1 Conceptual Questions 

The two-nucleon system is non-perturbative as evidenced by the presence 
of a shallow bound state (the deuteron) and large scattering lengths. Wein- 
berg [18 j showed that the strong enhancement of the scattering amplitude 
arises from purely nucleonic intermediate states. He therefore suggested to 
use perturbation theory to calculate the NN potential and to apply this 
potential in a scattering equation to obtain the NN amplitude. We adopt 
this prescription. 

Since the irreducible diagrams that make up the potential are calculated 



using covariant perturbation theory (cf. Sec. 5.1), it is consistent to start 
from the covariant Bethe-Salpeter (BS) equation [55] describing two-nucleon 
scattering. In operator notation, the BS equation reads 

T = V + VGT (81) 

with T the invariant amplitude for the two-nucleon scattering process, V the 
sum of all connected two-particle irreducible diagrams, and Q the relativistic 
two-nucleon propagator. The BS equation is equivalent to a set of two 
equations 

T = V + VgT (82) 
V = V + V{G-g)V (83) 
= V + V l7r {Q - g) V l7r + . . . , (84) 

where g is a covariant three-dimensional propagator which preserves rela- 
tivistic elastic unitarity. We choose the propagator g proposed by Blanken- 
becler and Sugar (BbS) jl3] (for more details on relativistic three-dimensional 



reductions of the BS equation, see Ref. [7]). The ellipsis in Eq. (84) stands 



for terms of irreducible 3tt and higher pion exchanges which we neglect. 
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Note that when we speak of covariance in conjunction with (heavy baryon) 
ChPT, we are not referring to manifest covariance. Relativity and relativis- 
tic off-shell effects are accounted for in terms of a Q/Mn expansion up to 
the given order. Thus, Eq. (1841) is evaluated in the following way, 



V « V(on-shell) + V ln Q V ln - V ln g V ln 



(85) 



where the pion-exchange content of V(on-shell) is V\ n + V 27T with Vi^ the on- 



shell 1PE given in Eq. ( 48 ) and V£ n the irreducible 2ir exchanges calculated 



in Sec. 5.1 but without the box. Vi n denotes the relativistic (off-shell) 1PE. 
Notice that the term (Vi^ Q Vi n — Vi n g Vi n ) represents what has been called 



"the (irreducible part of the) box diagram contribution" in Sec. 5.1 where 
it was evaluated at various orders. 

The full chiral NN potential V is given by irreducible pion exchanges 
V n and contact terms V c t, 

V = V n + (86) 



with 



Vlvr + V 2n + . . • , 



(87) 



where the ellipsis denotes irreducible 3ir and higher pion exchanges which 
are omitted. Two-pion exchange contributions appear in various orders 



V 



2tt 



V£> + V£> + V£> + .. 



'2tt 



r(3) 



2tt 



(4) 



2tt 



as calculated in Sec. 5.1. Contact terms come in even orders, 



V 



ct 



KS 0) + K ( t 2) + Ki 4) + 



(89) 



and were presented in Sec. 5.3 The potential V is calculated at a given order. 



(3) 

For example, the potential at NNLO includes 2PE up to V^ % and contacts 
up to ■ At N 3 LO, contributions up to V^ and v}^ are included. 



The potential V satisfies the relativistic BbS equation, Eq. (82 ). Defining 



V(p',p) 



1 M N M N 
V(p ,p) 



(2tt)3 V E p , 



E„ 



(90) 



and 



f(p',p) 



1 ¥» T f*> * Me. 



(91) 
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with E p = 



Mfj+p 2 (the factor l/(2vr) 3 is added for convenience 



the 



BbS equation collapses into the usual, nonrelativistic Lippmann-Schwinger 
(LS) equation, 



T(p',p) = V(p',p) + / d 3 p" V(p 



• i ->//\ 



M 



p" 2 + ie 



T(p",p). (92) 



Since V satisfies Eq. (92), it can be used like a usual nonrelativistic 



potential, and T is the conventional nonrelativistic T-matrix. 

Iteration of V in the LS equation requires cutting V off for high momenta 
to avoid infinities, This is consistent with the fact that ChPT is a low- 
momentum expansion which is valid only for momenta Q <C A x ~ 1 GeV. 
Thus, we multiply V with a regulator function 



V(p',p) 



V(p',p) e 



(p'/A) 2 « e -(p/A) 2 ' 1 



2n 



+ 



p\ 

A 



2n 



+ 



(93) 
(94) 



with the 'cutoff parameter' A around 0.5 GeV. Equation (94) provides an 



indication of the fact that the exponential cutoff does not necessarily affect 
the given order at which the calculation is conducted. For sufficiently large 
n, the regulator introduces contributions that are beyond the given order. 
Assuming a good rate of convergence of the chiral expansion, such orders are 
small as compared to the given order and, thus, do not affect the accuracy at 
the given order. In our calculations we use, of course, the full exponential, 



Eq. (93), and not the expansion. On a similar note, we also do not expand 



the square-root factors in Eqs. (90 91 ) because they are kinematical factors 



which guarantee relativistic elastic unitarity. 



5.4.2 What Order? 

Since in nuclear EFT we are dealing with a perturbative expansion, at some 
point, we have to raise the question, to what order of ChPT we have to go 
to obtain the precision we need. To discuss this issue on firm grounds, we 
show in Table [3] the x 2 /datum for the fit of the world np data below 290 
MeV for a family of np potentials at NLO and NNLO. The NLO potentials 
produce the very large x 2 /datum between 67 and 105, and the NNLO are 
between 12 and 27. The rate of improvement from one order to the other 
is very encouraging, but the quality of the reproduction of the np data at 
NLO and NNLO is obviously insufficient for reliable predictions. 
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Table 3: x 2 /datum for the reproduction of the 1999 np database [56] by 
families of np potentials at NLO and NNLO constructed by the Juelich 
group [57] . 



Tiab bin 


# of np 


- Juelich np potentials - 


(MeV) 


data 


NLO 


NNLO 


0-100 


1058 


4-5 


1.4-1.9 


100-190 


501 


77-121 


12-32 


190-290 


843 


140-220 


25-69 


0-290 


2402 


67-105 


12-27 



Based upon these facts, it has been pointed out in 2002 by Entem and 
Machleidt [251 I2S] that one has to proceed to N 3 LO. Consequently, the first 
N 3 LO potential was published in 2003 [27]. 

At N 3 LO, there are 24 contact terms (24 parameters) which contribute to 
the partial waves with L < 2 (cf. Sec.E3}. In Table H column 'Q 4 /N 3 LO', 



we show how these terms / parameters are distributed over the various partial 
waves. For comparison, we also show the number of parameters used in 
the Nijmegen partial wave analysis (PWA93) [49 and in the high-precision 
CD-Bonn potential [9]. The table reveals that, for S and P waves, the 
number of parameters used in high-precision phenomenology and in EFT at 
N 3 LO are about the same. Thus, the EFT approach provides retroactively 
a justification for what the phenomenologists of the 1990's were doing. At 
NLO and NNLO, the number of parameters is substantially smaller than for 
PWA93 and CD-Bonn, which explains why these orders are insufficient for 
a quantitative potential. This fact is also clearly reflected in Fig. [7j where 
phase shifts are shown for potentials constructed at NLO, NNLO, and N 3 LO. 



5.4.3 Charge-Dependence 

For an accurate fit of the low-energy pp and np data, charge-dependence 
is important. We include charge-dependence up to next-to-leading order of 
the isospin-violation scheme (NLO, in the notation of Ref. |58|). Thus, we 
include the pion mass difference in 1PE and the Coulomb potential in pp 
scattering, which takes care of the L0 contributions. At order NLO, we 
have the pion mass difference in 2PE at NLO, irj exchange [59|, and two 
charge-dependent contact interactions of order Q° which make possible an 
accurate fit of the three different 1 5'o scattering lengths, a pp , a nn , and a np . 
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Table 4: Number of parameters needed for fitting the np data in phase-shift 
analysis and by a high-precision NN potential versus the number of NN 
contact terms of EFT based potentials at different orders. 



Nijmegen CD-Bonn 
partial-wave high-precision 
analysis [49J potential [9] 


- Contact Potentials - 
Q° Q 2 Q 4 
LO NLO/NNLO N 3 LO 


% 3 4 
3 Si 3 4 


1 2 4 
1 2 4 


3 S 1 - 3 D 1 2 2 


1 3 


1 P 1 3 3 
3 P 3 2 
3 Pi 2 2 
3 P 2 3 3 


1 2 
1 2 
1 2 
1 2 


3 P 2 - 3 P 2 2 1 


1 


1 D 2 2 3 
3 I>1 2 1 
3 L> 2 2 2 
3 D 3 1 2 


1 
1 
1 
1 


3 L>3- 3 G , 3 1 





l F 3 1 1 
3 F 2 1 2 
3 F 3 1 2 
3 F 4 2 1 








3 F 4 - 3 ^4 





X G 4 1 
3 G 3 1 
3 G 4 1 
3 G 5 1 








Total 35 38 
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Figure 7: Phase parameters for np scattering as calculated from NN potentials at 
different orders of ChPT. The dotted line is NLO [57 , the dashed NNLO [57], and 
the solid N 3 LO [27] . Partial waves with total angular momentum J < 2 are dis- 
played. Solid dots represent the Nijmegen multienergy np phase shift analysis |49j 
and open circles are the GWU/VPI single-energy np analysis SM99 [50 . 



5.4.4 A Quantitative NN Potential at N 3 LO 

NN Scattering. The fitting procedure starts with the peripheral partial 
waves because they depend on fewer parameters. Partial waves with L > 3 
are exclusively determined by 1PE and 2PE because the N 3 LO contacts 
contribute to L < 2 only. 1PE and 2PE at N 3 LO depend on the axial- 
vector coupling constant, gA (we use qa = 1-29), the pion decay constant, 
/tt = 92.4 MeV, and eight low-energy constants (LECs) that appear in the 



dimension- two and dimension-three nN Lagrangians, Eqs. (37) and (38). 
In the fitting process, we varied three of them, namely, C2, C3, and C4. We 
found that the other LECs are not very effective in the NN system and, 
therefore, we kept them at the values determined from ttN (cf. Table [2]). 
The most influential constant is C3, which has to be chosen on the low side 
(slightly more than one standard deviation below its ttN determination) for 
an optimal fit of the NN data. As compared to a calculation that strictly 
uses the -irN values for C2 and C4, our choices for these two LECs lower 
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Table 5: x 2 /datum for the reproduction of the 1999 np database [56] 
by various np potentials. (Numbers in parentheses are the values of cutoff 
parameters in units of MeV used in the regulators of the chiral potentials.) 



Jlab bin 


# of np 


Idaho 


Juelich 


Argonne 


(MeV) 


data 


N 3 LO [27] 


N 3 LO [60] 


Vis m 






(500-600) 


(600/700-450/500) 




0-100 


1058 


1.0-1.1 


1.0-1.1 


0.95 


100-190 


501 


1.1-1.2 


1.3-1.8 


1.10 


190-290 


843 


1.2-1.4 


2.8-20.0 


1.11 


0-290 


2402 


1.1-1.3 


1.7-7.9 


1.04 


Table 6: x 2 /datum for the reproduction of the 1999 pp database [56] by 


various pp potentials. 


Notation as in Fig 






Ilab bin 


# of np 


Idaho 


Juelich 


Argonne 


(MeV) 


data 


N 3 LO [27] 


N 3 LO [60] 


^is m 






(500-600) 


(600/700-450/500) 




0-100 


795 


1.0-1.7 


1.0-3.8 


1.0 


100-190 


411 


1.5-1.9 


3.5-11.6 


1.3 


190-290 


851 


1.9-2.7 


4.3-44.4 


1.8 


0-290 


2057 


1.5-2.1 


2.9-22.3 


1.4 



the 3 i<2 and phase shifts bringing them into closer agreement with the 
phase shift analysis. The other F waves and the higher partial waves are 
essentially unaffected by our variations of c% and C4. Overall, the fit of all 
J > 3 waves is very good. 

We turn now to the lower partial waves. Here, the most important fit 
parameters are the ones associated with the 24 contact terms that contribute 
to the partial waves with L < 2. In addition, we have two charge-dependent 
contacts which are used to fit the three different 1 Sq scattering lengths, a pp , 
dnni and o> n p. 

In the optimization procedure, we fit first phase shifts, and then we 
refine the fit by minimizing the x 2 obtained from a direct comparison with 
the data. The x 2 /datum for the fit of the np data below 290 MeV is shown 
in Table [5] and the corresponding one for pp is given in Table [6] These tables 
show that at N 3 LO a \ 2 /datum comparable to the high- precision Argonne 
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Figure 8: Neutron-proton phase parameters as described by two potentials at 
N 3 LO. The solid curve is calculated from the Idaho N 3 LO potential [57] while the 
dashed curve is from the Juelich [50] one. Solid dots and open circles as in Fig. [7] 



Vis potential can, indeed, be achieved. The "Idaho" N 3 LO potential [27] 
produces a x 2 /datum = 1.1 for the world np data below 290 MeV which 
compares well with the x 2 /datum = 1.04 by the Argonne potential. In 2005, 
also the Juelich group produced several N 3 LO NN potentials [60 , the best 
of which fits the np data with a x 2 /datum = 1.7 and the worse with a 
X 2 /datum = 7.9 (see Table pi). While 7.9 is clearly unacceptable for any 
meaningful application, a x^datum of 1.7 is reasonable, although it does 
not meet the precision standard that few-nucleon physicists established in 
the 1990's. 

Turning to pp, the x 2 f° r PP data are typically larger than for np be- 
cause of the higher precision of pp data. Thus, the Argonne Vis produces a 
X 2 /datum = 1.4 for the world pp data below 290 MeV and the best Idaho 
N 3 LO pp potential obtains 1.5. The fit by the best Juelich N 3 LO pp poten- 
tial results in a x 2 /datum = 2.9 which, again, is not quite consistent with 
the precision standards established in the 1990's. The worst Juelich N 3 LO 
pp potential produces a x 2 /datum of 22.3 and is incompatible with reliable 
predictions. 
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Table 7: Deuteron properties as predicted by various NN potentials are com- 
pared to empirical information. (Deuteron binding energy Bd, asymptotic S state 
As, asymptotic D/S state r), deuteron radius rd, quadrupole moment Q, D-state 
probability Pd; the calculated r d and Q are without meson-exchange current con- 
tributions and relativistic corrections.) 





Idaho 


Juclich 










N 3 LO [27] 


N 3 LO [60] 


CD-Bonn [9J 


AV18[£U 


Empirical 




(500) 


(550/600) 








B d (MeV) 


2.224575 


2.218279 


2.224575 


2.224575 


2.224575(9) 


A s (fm- 1 / 2 ) 


0.8843 


0.8820 


0.8846 


0.8850 


0.8846(9) 


V 


0.0256 


0.0254 


0.0256 


0.0250 


0.0256(4) 


rd (fm) 


1.975 


1.977 


1.966 


1.967 


1.97535(85) 


Q (fm 2 ) 


0.275 


0.266 


0.270 


0.270 


0.2859(3) 


Pd (%) 


4.51 


3.28 


4.85 


5.76 




"See Table XVIII of Ref. [9] for references; the empirical value for r<j is 


from Ref. [62]. 



Phase shifts of np scattering from the best Idaho (solid line) and Juelich 
(dashed line) N 3 LO np potentials are shown in Figure |8j The phase shifts 
confirm what the corresponding \ 2 have already revealed. 

The Deuteron. The reproduction of the deuteron parameters is shown 
in Table [Tj We present results for two N 3 LO potentials, namely, Idaho [27] 
and Juelich [60] . Remarkable are the predictions by the chiral potentials for 
the deuteron radius which are in good agreement with the latest empirical 
value obtained by the isotope-shift method [62] ■ All AW potentials of the 
past (Table [7] includes two representative examples, namely, CD-Bonn [9] 
and AV18 [61]) fail to reproduce this very precise new value for the deuteron 
radius. 

In Fig. [9J we display the deuteron wave functions derived from the N 3 LO 
potentials and compare them with wave functions based upon conventional 
AW potentials from the recent past. Characteristic differences are notice- 
able; in particular, the chiral wave functions are shifted towards larger r 
which explains the larger deuteron radius. 

6 Many-Nucleon Forces 

As noted before, an important advantage of the EFT approach to nuclear 
forces is that it creates two- and many-nucleon forces on an equal footing. 
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Figure 9: Deuteron wave functions: the family of larger curves are S-waves, 
the smaller ones D- waves. The thick lines represent the wave functions 
derived from chiral NN potentials at order N 3 LO (thick solid: Idaho |27j . 
thick dashed: Juelich [60 ) . The thin dashed, dash-dotted, and dotted lines 
refer to the wave functions of the CD-Bonn [9], Nijm-I [8], and AV18 [61] 
potentials, respectively. 



6.1 Three-Nucleon Forces 



The first non-vanishing 3NF terms occur at NNLO and are shown in Fig. 10 



(cf. also Fig. [T] row 'Q 3 /NNLO', column '3N Force'). There are three di- 
agrams: the 2PE, 1PE, and 3N-contact interactions [39j HQ]. The 2PE 
3N-potential is given by 



T/3NF 
V 2PE 



9A 
2f. 



1 (&j ■ qi){aj ■ qj) 

tjg* (qf + ml)(q] + ml) 



F «/3 a B 
ijk ' t j 



(95) 



with qi = pi' — pi, where pi and p/ are the initial and final momenta of 
nucleon i, respectively, and 



7 a8 
ijk 



ra/3 



Aciml 2c 3 _ _ 



j - 



+ 



Qi ■ Qj 



f 2 
j - 



+ E ^ <?k ■ Wi x q 3 ] . (96) 



The vertex involved in this 3NF term is the two-derivative irirNN vertex 



(solid square in Fig. 10) which we encountered already in the 2PE contribu- 
tion to the iViV potential at NNLO. Thus, there are no new parameters and 
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Figure 10: The three-nucleon force at NNLO (from Ref. [30]). 



the contribution is fixed by the LECs used in NN. The 1PE contribution is 



« = D §- 2 £ • r,)^ • %) (97) 



rat 



and, finally, the 3N contact term reads 

V?r = EWrrr k . (98) 

The last two 3NF terms involve two new vertices (that do not appear in the 
2N problem) , namely, the ttNNNN vertex with parameter D and a 6 N ver- 
tex with parameters E. To pin them down, one needs two observables that 
involve at least three nucleons. In Ref. (ID], the triton binding energy and 
the nd doublet scattering length 2 a n d were used. Alternatively, one may also 
choose the binding energies of 3 H and 4 He [63] . Once D and E are fixed, the 
results for other 3N, 4N, . . .observables are predictions. In Refs. [Ml [63] . 
the first calculations of the structure of light nuclei ( 6 Li and 7 Li) were re- 
ported. Recently, the structure of nuclei with A = 10 — 13 nucleons has been 
calculated using the ab initio no-core shell model and applying chiral two 
and three-nucleon forces [65] . The results are very encouraging. Concerning 
the famous l A y puzzle', the above 3NF terms yield some improvement of 
the predicted nd A y , however, the problem is not solved [40J. 



Note that the 3NF expressions given in Eqs. (95 )-(98) above are the ones 
that occur at NNLO, and all calculations to date have included only those. 
Since we have to proceed to N 3 LO for sufficient accuracy of the 2NF, then 
consistency requires that we also consider the 3NF at N 3 LO. The 3NF at 



N 3 LO is very involved as can be seen from Fig. 11 but it does not depend 
on any new parameters. It is presently under construction [66]. So, for the 
moment, we can only hope that the A y puzzle may be solved by a complete 
calculation at N 3 LO. 
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Figure 11: Three- nucleon force contributions at N 3 LO (from Ref. [SB]). 



6.2 Four-Nucleon Forces 

In ChPT, four-nucleon forces (4NF) appear for the first time at N 3 LO 
(y = 4). Thus, N 3 LO is the leading order for 4NF. Assuming a good rate of 
convergence, a contribution of order (Q/A x ) 4 is expected to be rather small. 
Thus, ChPT predicts 4NF to be essentially insignificant, consistent with ex- 
perience. Still, nothing is fully proven in physics unless we have performed 
explicit calculations. Very recently, the first such calculation has been per- 
formed: The chiral 4NF, Fig. 12 has been applied in a calculation of the 
4 He binding energy and found to contribute a few 100 keV [HE]- It should 



be noted that this preliminary calculation involves many approximations, 
but it certainly provides the right order of magnitude of the result, which is 
indeed very small as compared to the full 4 He binding energy of 28.3 MeV. 



7 Conclusions 

The theory of nuclear forces has made great progress since the turn of the 
millennium. Nucleon-nucleon potentials have been developed that are based 
on proper theory (EFT for low-energy QCD) and are of high-precision, at 
the same time. Moreover, the theory generates two- and many-body forces 
on an equal footing and provides a theoretical explanation for the empirically 
known fact that 2NF > 3NF > 4NF .... 
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Figure 12: The four-nucleon force at N 3 LO (from Ref. [57]). 

At N 3 LO [26, 27J, the accuracy can be achieved that is necessary and 
sufficient for reliable microscopic nuclear structure predictions. First cal- 
culations applying the N 3 LO NN potential [27] in the conventional shell 
model [69j [70], the ab initio no-core shell model [7TJ [72] [73], the coupled 
cluster formalism [73] [75] [76] [77] [78] , and the unitary-model-operator ap- 
proach j79j have produced promising results. 

The 3NF at NNLO is known [39, 40^ and has been applied in few-nucleon 
reactions (30] [80] [81] as well as the structure of light nuclei [Ml E21 [65] . How- 
ever, the famous 'A y puzzle' of nucleon-deuteron scattering is not resolved 
by the 3NF at NNLO. Thus, one important outstanding issue is the 3NF at 
N 3 LO, which is under construction [66]. 

Another open question that needs to be settled is whether Weinberg 
power counting, which is applied in all current AW potentials, is consistent. 
This controversial issue is presently being debated in the literature [82] [83] . 
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A Fourth Order Two-Pion Exchange Contributions 



The fourth order 2PE contributions consist of two classes: the one-loop 
(Fig. [3]) and the two-loop diagrams (Fig. [4| . 

A.l One-loop diagrams 

This large pool of diagrams can be analyzed in a systematic way by intro- 
ducing the following well-defined subdivisions. 



A. 1.1 cf contributions. 

The only contribution of this kind comes from the football diagram with 
both vertices proportional to Cj (first row of Fig. One obtains |41j : 
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16tt 2 /^ 
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c\w 2 L(q) 
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(99) 



(100) 



A. 1.2 Cj/Mjv contributions. 



This class consists of diagrams with one vertex proportional to q and one 
l/M^r correction. A few graphs that are representative for this class are 
shown in the second row of Fig. [3] Symbols with a large solid dot and an 
open circle denote 1/Mjy corrections of vertices proportional to Cj. They are 
part of C n pf, Eq. (39). The result for this group of diagrams is 



Vc 

w c 

W T 



32ir 2 M N M 
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w LS 



c 4 L(g) 



.„ L(g) , 

^(Sm 2 + 5g 2 ) + w 2 



(104) 
(105) 



A. 1.3 1/Mjy corrections. 

These are relativistic 1/Mjy- corrections of the leading order 2-7T exchange 
diagrams. Typical examples for this large class are shown in row 3-6 of 
Fig. [3} T his time, there is no correction from the iterated 1PE, Eq. (65) or 
Eq. (66), since the expansion of the factor M^/E p does not create a term 
proportional to 1/Mjy. The total result for this class is 
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A. 2 Two-loop contributions. 

The two-loop contributions are quite intricate. In Fig. |4j we attempt a 
graphical representation of this class. The gray disk stands for all one-loop 
ttN graphs which are shown in some detail in the lower part of the figure. 
Not all of the numerous graphs are displayed. Some of the missing ones 
are obtained by permutation of the vertices along the nucleon line, others 
by inverting initial and final states. Vertices denoted by a small dot are 



from the leading order ttN Lagrangian C^L, Eq. (33), except for the four- 



pion vertices which are from C^J, Eq. (|27|). The solid square represents 



vertices proportional to the IjUjUs ck \v 

,(3) 



lich are introduced by the third 



ECs di 

order Lagrangian C)lIt, Eq. ( 381) . The d, vertices occur actually in one- 
loop iViV diagrams, but we list them among the two-loop NN contributions 
because they are needed to absorb divergences generated by one-loop ttN 
graphs. Using techniques from dispersion theory, Kaiser [H] calculated the 
imaginary parts of the NN amplitudes, Im V a (ifi) and Im W a (in), which 
result from analytic continuation to time-like momentum transfer q = ifi—0 + 
with [i > Im-x. From this, the momentum-space amplitudes V a {q) and 
W a {q) are obtained via the subtracted dispersion relations: 



v c ,s(q) 

V T (q) 



2q e 



dp, 



Im V c ,s(i^) 



vr hm* /^ 5 (^ 2 + q 2 ) ' 



2q 4 



dfx 



Im Vr(i[J>) 



vr J2m^ fi 3 (n 2 + q 2 ) ' 



(113) 
(114) 



and similarly for Wc,s,T- 

In most cases, the dispersion integrals can be solved analytically and the 
following expressions are obtained [26 : 
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+ 
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where k = y/ // 2 /4 — m 2 . 

Note that the analytic solutions hold modulo polynomials. We have 
checked the importance of those contributions where we could not find an 
analytic solution and where, therefore, the integrations have to be performed 
numerically. It turns out that the combined effect on NN phase shifts from 
^,l#\andvf 
than 0.01 deg in H waves, at T[ a b 

This renders these contributions negligible. Therefore, we omit W§\ V± b) , 
and Vg in the construction of chiral NN potentials at order N 3 LO. 



s is smaller than 0.1 deg in F and G waves and smaller 



300 MeV (and less at lower energies). 



In Eqs. (117) and (121), we use the scale-independent LECs, di, which 
are obtained by combining the scale-dependent ones, d T -(X), with the chiral 
logarithm, ln(m^/A), or equivalently di = dl(m n ). The scale-dependent 
LECs, (i[(A), are a consequence of renormalization. For more details about 
this issue, see Ref. [37] . 



B Partial Wave Decomposition of the Fourth Or- 
der Contact Potential 



The contact potential contribution of order four, Eq. (80), decomposes into 
partial- waves as follows. 
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The coefficients in the above expressions are given by: 
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-^D 9 + ^D 10 -^D 1B (126) 

References 

[1] H. Yukawa, Proc. Phys. Math. Soc. Japan 17, 48 (1935). 

[2] Prog. Theor. Phys. (Kyoto), Supplement 3 (1956). 

[3] M. Taketani, S. Machida, and S. Onuma, Prog. Theor. Phys. (Kyoto) 
7, 45 (1952). 

[4] K. A. Brueckner and K. M. Watson, Phys. Rev. 90, 699; 92, 1023 
(1953). 

[5] A. R. Erwin et al, Phys. Rev. Lett. 6, 628 (1961); B. C. Maglic et al, 
ibid. 7, 178 (1961). 

[6] Prog. Theor. Phys. (Kyoto), Supplement 39 (1967); R. A. Bryan and 
B. L. Scott, Phys. Rev. 177, 1435 (1969); M. M. Nagels et al, Phys. 
Rev. D 17, 768 (1978). 

[7] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[8] V. G. J. Stoks et al, Phys. Rev. C 49, 2950 (1994). 

[9] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

[10] A. D. Jackson, D. O. Riska, and B. Verwest, Nucl. Phys. A249, 397 
(1975). 

[11] R. Vinh Mau, in Mesons in Nuclei, edited by M. Rho and D. H. Wilkin- 
son (North- Holland, Amsterdam, 1979), Vol. I, p. 151. 

[12] M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires, 
and R. de Tourreil, Phys. Rev. C 21, 861 (1980). 

[13] R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

[14] F. Myhrer and J. Wroldsen, Rev. Mod. Phys. 60, 629 (1988). 

[15] D. R. Entem, F. Fernandez, and A. Valcarce, Phys. Rev. C 62, 034002 
(2000). 



50 



[16] G. H. Wu, J. L. Ping, L. J. Teng, F. Wang, and T. Goldman, Nucl. 
Phys. A673, 273 (2000). 



[17; 

[18 

[19 

[20 
[21 

[22 
[23 

[24 
[25 
[26 

[27; 

[28 

[29 

[30 
[31 
[32 

[33 
[34 



S. Weinberg, Physica 96A, 327 (1979). 

S. Weinberg, Phys. Lett. B 251, 288 (1990); Nucl. Phys. B363, 3 
(1991); Phys. Lett. B 295, 114 (1992). 

C. Ordohez, L. Ray, and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994); 
Phys. Rev. C 53, 2086 (1996). 

U. van Kolck, Prog. Part. Nucl. Phys. 43, 337 (1999). 

L. S. Celenza et al, Phys. Rev. C 46, 2213 (1992); C. A. da Rocha et 
al, ibid. 49, 1818 (1994); D. B. Kaplan et al, Nucl. Phys. B478, 629 
(1996). 

N. Kaiser, R. Brockmann, and W. Weise, Nucl. Phys. A625, 758 (1997). 

N. Kaiser, S. Gerstendorfer, and W. Weise, Nucl. Phys. A637, 395 
(1998). 

E. Epelbaum et al, Nucl. Phys. A637, 107 (1998); A671, 295 (2000). 

D. R. Entem and R. Machleidt, Phys. Lett. B 524, 93 (2002). 

D. R. Entem and R. Machleidt, Phys. Rev. C 66, 014002 (2002). 

D. R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 (2003). 

R. Machleidt and D. R. Entem, J. Phys. G: Nucl. Phys. 31, S1235 
(2005). 

P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 52, 339 
(2002). 

S. Scherer and M. R. Schindler, arXiv:hep-ph/0505265. 

Review of Particle Physics, J. Phys. G: Nucl. Part. Phys. 33, 1 (2006). 

S. Coleman, J. Wess, and B. Zumino, Phys. Rev. 177, 2239 (1969); C. 
G. Callan, S. Coleman, J. Wess, and B. Zumino, ibid. 177, 2247 (1969). 

J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984). 

J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B307, 779 (1988). 



51 



[35] V. Bernard, N. Kaiser, and U.-G. Meifiner, Int. J. Mod. Phys. E 4, 193 
(1995). 

[36] N. Fettes, U.-G. Meifiner, M. Mojzis, and S. Steininger, Ann. Phys. 
(N.Y.) 283, 273 (2000); 288, 249 (2001). 

[37] N. Fettes, U.-G. Meifiner, and S. Steiniger, Nucl. Phys. A640, 199 
(1998). 

[38] N. Kaiser, Phys. Rev. C 61, 014003 (1999); 62, 024001 (2000). 

[39] U. van Kolck, Phys. Rev. C 49, 2932 (1994). 

[40] E. Epelbaum et al, Phys. Rev. C 66, 064001 (2002). 

[41] N. Kaiser, Phys. Rev. C 64, 057001 (2001). 

[42] N. Kaiser, Phys. Rev. C 65, 017001 (2002). 

[43] R. Blankenbecler and R. Sugar, Phys. Rev. 142, 1051 (1966). 

[44] This section closely follows Ref. [26] . 

[45] G. Q. Li and R. Machleidt, Phys. Rev. C 58, 3153 (1998). 

[46] V. Stoks, R. Timmermans, and J. J. de Swart, Phys. Rev. C 47, 512 

(1993) . 

[47] R. A. Arndt, R. L. Workman, and M. M. Pavan, Phys. Rev. C 49, 2729 

(1994) . 

[48] P. Biittiker and U.-G. Meifiner, Nucl. Phys. A668, 97 (2000). 

[49] V. G. J. Stoks, R. A. M. Klomp, M. C. M. Rentmeester, and J. J. de 
Swart, Phys. Rev. C 48, 792 (1993). 

[50] R. A. Arndt, I. I. Strakovsky, and R. L. Workman, SAID, Scattering 
Analysis Interactive Dial-in computer facility, George Washington Uni- 
versity (formerly Virginia Polytechnic Institute), solution SM99 (Sum- 
mer 1999); for more information see, e. g., R. A. Arndt, I. I. Strakovsky, 
and R. L. Workman, Phys. Rev. C 50, 2731 (1994). 

[51] In fact, preliminary calculations, which take an important class of dia- 
grams of order five into account, indicate that the N 4 LO contribution 
may prevailingly be repulsive (N. Kaiser, private communication). 



52 



[52] G. E. Brown and A. D. Jackson, The Nucleon-Nucleon Interaction, 
(North-Holland, Amsterdam, 1976). 

[53] N. Kaiser, Phys. Rev. C 63, 044010 (2001). 

[54] K. Erkelenz, R. Alzetta, and K. Holinde, Nucl. Phys. A176, 413 (1971); 
note that there is an error in equation (4.22) of this paper where it 
should read 



2J-1 

and 



+ WL = 2qq' J + 2 r^+2.(Q)_ a AO) 
Lb HH 2J + 3 



[55 
[56 

[57; 

[58 
[59 
[60 

[61 
[62 
[63 

[64 
[65 

[66 

[67; 



E. E. Salpeter and H. A. Bethe, Phys. Rev. 84, 1232 (1951). 
The 1999 NN data base is defined in Ref. [9]. 

E. Epelbaum, W. Glockle, and U.-G. Meifiner, Eur. Phys. J. A19, 401 

(2004) . 

M. Walzl et al, Nucl. Phys. A693, 663 (2001). 

U. van Kolck et al, Phys. Rev. Lett. 80, 4386 (1998). 

E. Epelbaum, W. Glockle, and U.-G. Meifiner, Nucl. Phys. A747, 362 

(2005) . 

R. B. Wiringa et al, Phys. Rev. C 51, 38 (1995). 
A. Huber et al, Phys. Rev. Lett. 80, 468 (1998). 

A. Nogga, P. Navratil, B. R. Barrett, and J. P. Vary, Phys. Rev. C 73, 
064002 (2006). 

A. Nogga et al, Nucl. Phys. A737, 236 (2004). 

P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand, and A. 
Nogga, |arXiv:nucl-th/0701038[ 

U.-G. Meifiner, Proc. 18th International Conference on Few-Body Prob- 
lems in Physics, Santos, SP, Brazil, August 2006, to be published in 
Nucl. Phys. A. 

E. Epelbaum, Phys. Lett. B 639, 456 (2006). 



53 



[68] D. Rozpedz ik et al, Acta Phys. Polon. B37, 2889 (2006); |arXiv:nucl- 
|th/06060l7l 



[69 
[70 
[71 
[72 
[73 
[74 
[75 
[76 

[77; 
[78 
[79 
[80 
[81 

[82 
[83 



L. Coraggio et al, Phys. Rev. C 66. 021303 (2002). 
L. Coraggio et al, Phys. Rev. C 71. 014307 (2005). 
P. Navratil and E. Caurier (2004) Phys. Rev. C 69 014311. 

C. Forssen et al, Phys. Rev. C 71, 044312 (2005). 
J.P. Vary et al, Eur. Phys. J. A 25 sOl, 475 (2005). 
K. Kowalski et al, Phys. Rev. Lett. 92, 132501 (2004). 

D. J. Dean and M. Hjorth-Jensen (2004) Phys. Rev. C 69 054320. 

M. Wloch et al, J. Phys. G 31, S1291 (2005); Phys. Rev. Lett. 94, 
21250 (2005). 

D. J. Dean et al, Nucl. Phys. 752, 299 (2005). 
J.R. Gour et al, Phys. Rev. C 74, 024310 (2006). 

S. Fujii, R. Okamato, and K. Suzuki, Phys. Rev. C 69, 034328 (2004). 

K. Ermisch et al, Phys. Rev. C 71, 064004 (2005). 

H. Witala, J. Golak, R. Skibinski, W. Glockle, A. Nogga, E. Epelbaum, 
H. Kamada, A. Kievsky, and M. Viviani, Phys. Rev. C 73, 044004 
(2006). 

A. Nogga, R. G. E. Timmermans, and U. van Kolck, Phys. Rev. C 72, 
054006 (2005). 

E. Epelbaum, U.-G. Meifiner, |a7Xiv:nucl-th/ 0609037, 



54 



